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Random walk models, such as the trap model, continuous time random walks, and comb models 
exhibit weak ergodicity breaking, when the average waiting time is infinite. The open question is: 
what statistical mechanical theory replaces the canonical Boltzmann-Gibbs theory for such systems? 
In this manuscript a non-ergodic equilibrium concept is investigated, for a continuous time random 
walk model in a potential field. In particular we show that in the non-ergodic phase the distribution 
of the occupation time of the particle on a given lattice point, approaches U or W shaped distribu- 
tions related to the arcsin law. We show that when conditions of detailed balance are applied, these 
distributions depend on the partition function of the problem, thus establishing a relation between 
the non-ergodic dynamics and canonical statistical mechanics. In the ergodic phase the distribution 
function of the occupation times approaches a delta function centered on the value predicted based 
on standard Boltzmann-Gibbs statistics. Relation of our work with single molecule experiments is 
briefly discussed. 

PACS numbers: 05.20.-y, 05.40.-a, 02.50.-r, 05.90.+m 



I. INTRODUCTION 

There is growing interest in non-ergodicity of systems 
whose dynamics is governed by power law waiting times, 
is such a way that a state of the system is occupied 
with a sojourn time whose average is infinite. Such non- 
ergodicity, called weak ergodicity breaking 0], was first 
introduced in the context of glassy dynamics. It has 
found several applications in Physics: phenomenological 
models of glassy dynamics pj, laser cooling 0, blinking 
quantum dots ]3J, , and models of atomic transport in 
optical lattice [5|- For example single blinking quantum 
dots, when interacting with a continuous wave laser field, 
turn at random times from a bright state in which many 
fluorescent photons are emitted, to a dark state. It is 
found that the distribution of dark and bright times fol- 
lows power law behavior. Somewhat similar statistical 
behavior is found also for laser cooling of atoms, where 
the atom is found in two states in momentum space, a 
cold trapped state and a free state, the sojourn time 
probability density function has a power law behavior 
ip(r) oc T~( 1+a } with a < 1. For such systems the time 
average of physical observable, for example the time av- 
erage of fluorescence intensity of single quantum dots, is 
non-identical to the ensemble average even in the long 
time limit. From a stochastic point of view such ergod- 
icity breaking is expected, since the condition to obtain 
ergodicity is that the measurement time t be much longer 
than the microscopical time scale of the problem. How- 
ever the microscopic time scale in our examples is infi- 
nite, namely the mean trapping times or the mean dark 
and bright times diverge. When these characteristic time 
scales are infinite, namely a < 1, we can never make time 
averages for long enough times to obtain ergodicity. 

It is important to note that the concept of a waiting 
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(i.e. trapping) time probability density function (PDF) 
^(t), with diverging first moment, is wide spread and 
found in many fields of Physics 0, S S H |E E3- It 
was introduced into the theory of transport of charge 
carriers in disordered material in the context of the 
continuous time random walk (CTRW). The CTRW de- 
scribes a random walk on a lattice with a waiting time 
PDF of times between jump events i/j(t). The model 
exhibits anomalous diffusion and aging behaviors 
[lH IrH 111 . [Tg| . when a < 1, which are related to er- 
godicity breaking. The CTRW found many applications 
in the context of chaotic dynamics 0,0,0, tracer dif- 
fusion in complex flows 0, 0] , financial time series [2(| , 
diffusion of bead in a polymer network |2l| . to name a 
few H 0. The dynamics of CTRW is similar to the 
dynamics of the comb model [8( and the annealed ver- 
sion of the trap model |16(. In turn the trap model is 
related to the random energy model [23J. All these sys- 
tems and models can be at-least suspected of exhibiting 
non-ergodic behavior, and hence constructing a general 
theory of non-ergodicity for such systems is in our opin- 
ion a worthy goal. 

Systems and models exhibiting anomalous diffusion, 
and CTRW behaviors can be divided into two categories. 
Systems where the random walk is close to thermal equi- 
librium, where the temperature of the system is well de- 
fined at least from an experimental point of view, and 
non-thermal systems. The ergodicity breaking of ther- 
mal CTRW models, is in conflict with Boltzmann-Gibbs 
ergodic assumption. As far as we know, there is no theory 
characterizing the non-ergodic properties of the CTRW 
for either thermal or non-thermal type of random walk. 

Hence one goal of this manuscript is to obtain the non- 
ergodic properties of the well known CTRW model on a 
lattice. Secondly we investigate ergodicity breaking and 
its relation to Boltzmann-Gibbs statistics. Using rather 
general arguments and using a CTRW model we inves- 
tigate the distribution of the total occupation times of 
a lattice point or a state of the system. We show that 
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in the limit of long measurement time and in the er- 
godic phase the occupation times are obtained using the 
Boltzmann-Gibbs canonical ensemble, provided that de- 
tailed balance conditions are satisfied. In the non-ergodic 
phase we obtain non-trivial distribution of the occupa- 
tion times, which are related to the arcsin law. These 
limiting distributions are unique in the sense that they 
do not depend on all the dynamical details of the under- 
lying model. Further the distributions we obtain depend 
on Boltzmann's probability namely on the temperature 
T and the partition function Z . Thus a connection is 
established between non-ergodic dynamics and the basic 
tool of statistical mechanics. 

The study of occupation times in the context of clas- 
sical Brownian motion was considered by P. Levy. Con- 
sider a Brownian path generated with x(t) = rj(t), where 
r)(t) is Gaussian white noise, in the time interval (0,i), 
and with free boundary conditions. The total time t+, 
the particle spend on the half space x > is called the 
occupation time of the positive half space. The fraction 
of occupation time p + = t +/t is distributed according to 
the celebrated arcsin law 123 



lim / (p+) 



1 



TT^/p+ (1 -p+) 



(1) 



where < p + < 1. In contrast to naive expectation, 
it is unlikely to find p + — 1/2, which would mean that 
the particle remains half of the time in i > 0. Instead 
f (p + ) diverges on p + = and p~ = 1, indicating that 
the Brownian particle tends to stay either in x > or 
in x < for long times of the order of the measurement 
time t. Hence f{p + ) has a U shape. Such a behavior is 
related to the survival probability of the Brownian par- 
ticle. The probability of a Brownian particle, starting at 
x > to remain in x > without crossing x = 0, decays 
like a power law t- 1 ' 2 . The average time the particle 
remains in i > 0, before the first crossing of x = 0, is 
infinite. Similar U shape distributions, in far less trivial 
examples, are investigated more recently in the context 
of random walks in random environments |25j . renewal 
processes [2(j , stochastic processes zero temperature 
Glauber spin dynamics 1281. diffusion equation |31| . two 
dimensional Ising model p30| , and growing interface |29| . 

The study of non-ergodicity within the CTRW frame- 
work is timely due to recent single molecule ]^ type 
of experiments. In many experiments anomalous diffu- 
sion, and power law behavior was observed using single 
particle tracking techniques 0, El 03 H3 ( e -S- sm gl e 
quantum dots |4|). An interesting example is the diffusive 
motion of magnetic beads in an actin network |2l|. The 
latter exhibit a CTRW type of behavior while the system 
has a well defined temperature T, namely the random 
walk seems close to thermal equilibrium and the particle 
is coupled to a thermal heat bath. In particular, long 
tailed waiting time distributions were recorded 

and anomalous sub-diffusion (r 2 ) ~ t a with a < 1, was 
observed. While clearly ensemble average classification 
of the anomalous process, e.g. the mean square displace- 



ment, are important, it is the time averages of single 
particle trajectories which distinguish the single parti- 
cle measurement from standard ensemble average type of 
measurement. And stochastic theories of non-ergodicity 
can help with the fundamental question in single molecule 
experiments: are time averages recorded in such experi- 
ment identical to the corresponding ensemble averages? 
and if not how do we classify the non-ergodic phase? 

This paper is organized as follows. In Sec. [H]we discuss 
a possible generalization of Boltzmann-Gibbs statistics 
for non-ergodic dynamics. In Sec. IIIII we introduce the 
CTRW model, which yields the non-ergodic dynamics. 
Sec. II VI is the main technical part of the paper, in which 
we obtain first passage time properties of the CTRW. 
The relation of these properties to the non-ergodic be- 
havior is shown. In Sec. [V]we give the main results and 
compare between the non-ergodic framework and stan- 
dard Boltzmann-Gibbs statistics. A brief summary of 
our results was published recently ■ 



II. FROM BOLTZMANN STATISTICS TO 
NON-ERGODICITY 

In this section we discuss a possible non-ergodic gener- 
alization of Boltzmann-Gibbs theory, without attempt- 
ing to prove its validity. 

The basic tool in statistical mechanics is Boltzmann's 
probability P x of finding a system in a state with energy 
E x , 



P B = 



exp (-§f) 



Z 



(2) 



where T is the temperature and Z — ^2 x exp(—E x /T). 
In Eq. J2J we use the canonical ensemble and assume a 
classical system, with discrete energy states {0 < E\ < 
E% • • •}. To obtain the average energy of the system, we 
use 



(e) = J2e x p^ 



(3) 



and similarly for other physical observables like entropy, 
free energy etc. Eq. J2J) is an ensemble average. When 
measurement of a single system is made, a time aver- 
age of a physical observable is recorded. Consider a sys- 
tem randomly changing between its energy states {E x }. 
At a given time the system occupies one energy state. 
Let t x be the total time spent by the system in energy 
state E x , within the total observation period (0, i). The 
system may visit state E x many times during the evolu- 
tion, hence t x is composed in principle from many sojourn 
times. We define the occupation fraction 

p x = —, and the time average energy is E = E x p x . 

X 

(4) 

According to statistical mechanics, once the ergodic hy- 
pothesis is satisfied, and within the canonical formalism 
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Px = Px anc ^ then E = (E) , and similarly for other phys- 
ical observables. More generally, the occupation fraction 
p x is a random variable, whose statistical properties de- 
pend on the underlying dynamics. If Boltzmann's condi- 
tions hold the probability density function (PDF) of p x 



/(Px) 



[Px 



P l 



(5) 



in the thermodynamic limit. The last Eq. is a restate- 
ment of the ergodic hypothesis. 

In this manuscript we discuss a possible generalization 
of the ergodic hypothesis. Our proposal is that the PDF 
of p x , for certain models described by CTRW type of 
dynamics, is described by a S a function 

f(Vx) = S a (K x ,p x ) = 



KxPT 1 (1 



_ \a~ 
Px) 



* K x 2 (1- p x ) 2a +p 2 x a + 2U X (1- p x fp a x cosTra ' 

(6) 

This PDF was obtained by Lamperti 36] in the context 
of the mathematical theory of occupation times (and see 
Appendix A for details). For 1Z X = I, a = 1/2 we have 
the arcsin law. Here we claim that when local detailed 
balance condition is satisfied 



n x 



i-p?'- 



(7) 



and < a < 1. When a=lwe get usual ergodic be- 
havior defined in Eq. JSJ. Eq. ijfJJl is valid only in the 
limit of long measurement time. In the non-ergodic phase 
a < 1 Eqs. (101711 establish a relation between the er- 
godicity breaking and Boltzmann-Gibbs statistics. The 
exponent a is the anomalous diffusion exponent in the 
relation (x 2 ) oc t a . 

For CTRWs not satisfying detailed balance condition 
a more general rule holds. We will show that the PDF of 
the fraction of time spent on lattice point x, p x , is still 
given by Eq. ©. However now 



n, 



peq 



l-P 



eq • 



(8) 



where P x q is the probability that a particle occupies lat- 
tice point x in equilibrium (an equilibrium is obtained for 
system of finite size). Here P x q and P x are probabilities 
in ensemble sense, namely if we consider an ensemble of 
N non interacting particles (or systems) satisfying some 
dynamical rule, P x q and P x yield in principle the proba- 
bility that a member of the ensemble occupies state x in 
equilibrium, which is not identical to p x for non-ergodic 
systems. 

Let us give some general arguments for the validity of 
Eqs. (|6I7|I . Consider a particular energy state of the 
system and call it E x . At a given time the system is 
either in energy state E x or is in any of the other energy 
states. When the system does not occupy state x we will 



say that the system is in state nx (not x). Assume that 
sojourn times in states x and nx are 



^x{t) 



\T{-a)\r 



l+a 



tpnx(T) 



A,, 



-a t 



i+ 



a ( 9 ) 



when r is large. Also assume that sojourn times in states 
x and nx are not correlated. Thus we imagine the system 
occupying state x then occupying state nx, then again 
state x etc. The amplitudes A x and A nx will generally de- 
pend on the particular dynamics of the system. We show 
in Appendix A that Eq. © holds with 1Z X = A x /A nx . 
Generally it seems a hopeless mission to calculate the 
ratio A x /A nx from any microscopical model. However 
a simple physical argument yields the ratio 1Z X . Assume 
that for an ensemble of systems Boltzmann-Gibbs statis- 
tical mechanics holds. Such an assumption means that 
on average we must have 



(P, 



Pi 



(10) 



where P x is Boltzmann's probability of finding a member 
of an ensemble of systems in state x. On the other hand, 
Eq. © yields 



(Px 



Pxf(Px) d Px 



(ii) 



Using Eq. jTUI EU we obtain Eq. 0. 

Our work, is related to the concept of weak ergodicity 
breaking, suggested by Bouchaud 1]. In standard statis- 
tical mechanics, one divides the phase space of the sys- 
tem, into equally sized cells, and the system is supposed 
to visit these cells, with equal probability under certain 
constrains (e.g. the energy of the system is constant for 
the micro-canonical ensemble). Strong ergodicity break- 
ing means that in order to leave one phase space cell 
to another, one has to cross a barrier (e.g. an energy 
barrier) which becomes infinite, in the thermodynamic 
limit. In this case the time it takes for the system to 
move from one state to the other is infinite. It is worth 
while thinking of such a process, in terms of a distribu- 
tion of escape times, ip(r) = i?exp(— Rt), where Rt is 
small and t is the measurement time (e.g. an activation 
over a very high energy barrier). In that case the parti- 
cle/system simply remains in a certain domain of phase 
space, for the whole period of observation, and the sys- 
tem does not explore its entire phase space available for 
ergodic systems. A very different scenario was suggested 
by Bouchaud, in the context of glassy dynamics and the 
trap model. If the distribution of sticking times, follows 
power law behavior, the average escape time diverges 



- a dT = 



(12) 



when a < 1. Note that also for the strong non-ergodicity 
case we may have an infinite waiting time (r) = 1/R 
when R — > 0. However for power law waiting times 
the system or particle may still explore its phase space. 
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Or in other words, exponential waiting times and power 
law waiting times, yield very different type of dynam- 
ics, even if for both the average waiting time is infinite. 
Thus, roughly speaking, for weak non-ergodicity and for 
ensemble of particles we may still get Boltzmann-Gibbs 
statistics, since from any initial condition the phase space 
is totally covered. However the system remains weakly 
non-ergodic, since during its evolution, the system will 
randomly pick one state, which it will occupy for a very 
long period (but it still visits all the other states) and 
then time averages are not equal to ensemble averages. 
The goal of this manuscript is to show that the strong 
assumptions we used are correct within a specific model, 
the well known CTRW model. 



III. CTRW IN FORCE FIELD 

We consider a one dimensional CTRW walk on a lat- 
tice. The lattice points are labeled with index x and 
x = —L,—L + 1, 0, ...L, hence the system size is 
2L + 1. On each lattice point we define a probability 
< Qr{x) < 1 for jumping right, and a probability for 
jumping left Ql{x) = 1 — Qr(x). Let ip(r) be the PDF 
of waiting times at the sites, this PDF does not depend 
on the position of the particle. If the particle starts at 
site x = 0, it will wait there for a period t\ determined 
from iP(t), it will then jump with probability Ql(0) to 
the left, and with probability Qr(0) to the right. After 
the jump, say to lattice point 1, the particle will pause for 
a period T2, whose statistical properties are determined 
by ip(r). It will then jump either back to point or to 
x = 2, according to the probability law Qr(1). Then 
the process is renewed. We consider reflecting boundary 
conditions, namely Ql{L) = Qr(—L) = 1. 

The case of a long tailed waiting time distribution, 
where ip(r) cx r~( 1+a ' > when r — ► oo and < a. < 1 
yields a non-ergodic behavior. In this case the average 
waiting time is infinite. The Laplace transform of ip(r) 
is 



i> (u) = / e~ UT ip (t) dr. 



o 

As-usual according to Tauberian theorem ^4 
u behavior is 

ip(u) ~ 1 - Au a + ■■■ 



(13) 



the small 



(14) 



and A > is a constant. 

Choose a specific lattice point x, then define 9 x (t) = 1 
if the particle is on x, otherwise it is zero. We define the 
occupation fraction as the time average of 9 x (t), 



P.r 



t 



(15) 



namely p x — t x /t where t x is the total time spent on 
lattice point x (i.e., the occupation time of site x ). We 
will later calculate the PDF of p x . 



Two special cases are the unbiased CTRW, where 
Q L (x) = Q R {x) = 1/2, and the biased CTRW with 
Ql(x) = q. In these cases all transition probabilities 
do not depend on the position of the random walker x, 
besides on the boundaries of-course. In the language of 
random walks these cases describe symmetric diffusion 
process, and diffusion with a drift. Note that in our 
model Ql(x) are not random variables, rather they are 
included in the model to mimic a deterministic poten- 
tial field acting on the system. For detailed discussion of 
CTRW models see [10. 

The case of diffusion with a constant drift, i.e., q ^ 1/2 
is used many times to model diffusion under the influence 
of a constant external driving force T . If the Physical 
process is close to thermal equilibrium the condition of 
detailed balance is imposed on the dynamics, in order 
that for an ensemble of particles Boltzmann equilibrium 
is reached [see further discussion after Eq. i|18fl ]. The 
potential energy at each point x, due to the interaction 
with the external driving force is E(x) — —Tax and a 
is the lattice spacing. The condition of detailed balance 
then reads 



Qr{x) 



exp 



J- a 



(16) 



where T is the temperature, and the right hand side 
of Eq. I|16fl is independent of lattice coordinate. Since 
Ql{x) = q is independent of x we have 



1 



q = 



1 + cxp(Fa/T) 



(17) 



More generally we define an energy profile for the sys- 
tem {E-l, E-l+i, —,Ei, ■ ■ ■}■ The general detailed bal- 
ance condition is then 



l-Q L (x-l) 



exp 



E x -i — E x 
T 



(18) 



The choice of detailed balance condition means that 
for an ensemble of particles standard Boltzmann-Gibbs 
statistics holds. Thus for example if we observe many 
independent particles, and look at their density profile in 
equilibrium, we will see a profile which is determined by 
Boltzmann equilibrium. On the other hand if we consider 
a trajectory of a single particle, and from it find p x we 
are not likely to find the value of p x close to Boltzmann's 
probability, when a < 1. Thus ergodicity breaking is 
found on the level of a single particle. Note that there 
is an interesting transition between one particle informa- 
tion and many particle behavior, however this is not the 
subject of our work |31| - 



IV. FIRST PASSAGE TIMES 

The problem of ergodicity breaking is related in this 
section to the problem of first passage times. 
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The process 9 x (t) is a two state process, with state x 
denoting particle on lattice point x and state nx indicat- 
ing that the particle is not on x. Obviously the waiting 
times in state x are given by ip x (r) = ip( T )- To obtain the 
PDF of waiting times in state nx, ip n w( T ) we must cal- 
culate statistical properties of first passage times. After 
the particle leaves point x it is located either on x + 1 or 
x — 1 with probabilities Qr(x) and Ql(x), respectively. 
Let tL denote the time it will take the particle to return 
to x starting at point x — 1, i.e. the first passage time 
from x — 1 to x. Let t R be the first passage time to reach 
x starting from x + 1. Let fR(t R ) be the PDF 

of the first passage time tu \ti] respectively. Then the 
PDF of times in state nx is given by 



V>n* (t) = Q R (x)f R (r) + Q L {x)f L (r) 



In principle once the long time behavior of the PDFs 
of first passage times is obtained, we have tp nx (r) and 
^ x (r), and then we may use the formalism developed in 
Appendix A to obtain the PDF of the occupation frac- 
tion p x . We now investigate the first passage times PDFs 
for biased and unbiased CTRW, using an analytical ap- 
proach. The reader not interested in mathematical de- 
tails, may skip to Sec. |Vj 



A. Relation Between Discrete Time and 
continuous time RWs 

For convenience we define a new lattice. We con- 
sider the CTRW in one dimension, on lattice points 
x = 0, 1,2, ...,JD. Point x = is a "sticky" absorbing 
boundary, namely once the particle reaches point x = 
it remains there for ever. Point L is a reflecting bound- 
ary, and initially at time t = the particle is on x = 1. 
Let Sqt (t) be the survival probability of the CTRW par- 
ticle, and the subscript CT indicates CTRW. The object 
of interest is the PDF of first passage time Jct (t) , which 
is minus the time derivative of <Sct(£)- The solution is 
possible due to an important relation between the 
CTRW first passage time problem and that of discrete 
time random walks. In [37j first passage time problem 
with CTRW dynamics with exponential waiting times 
was considered. 

Point of the new lattice is point x in the original 
problem and L = L — x and similarly for the other L — 1 
points of the new lattice. Hence the calculation of the 
first passage PDF on the new lattice x = 0, 1, 2, ...L yields 
faita)- With straight forward change of notation we may 
consider also 

Let S'ci'(i) be the survival probability of the CTRW 
particle in the interval x = 1, ....,x = L. Let S^ S (N) 
be the probability of survival after N jumps events, for 
a particle starting at x — 1, the subscript dis stands for 
discrete. Then 



(20) 



where P(N, t) is the probability for N steps, in time t, in 
a CTRW process. In Laplace, t — *■ u space it is easy to 
show using the convolution theorem of Laplace transform 
that 



P(N,u) = 



1 - tM in/ n 



(21) 



where ip(u) is the Laplace transform of ^P{t). In this work 
the discrete Laplace transform of an arbitrary function 
G(N), also called the z transform is defined as 



G(z)= J2z N G(N) 



N=0 



(19) Using Eqs. QXWfy we find 



1 - i>{u) ~ 



^dis 



(22) 



(23) 



This equation establishes the relation between the dis- 
crete and continuous time problems. 

Let P x (N) be the probability of occupying site x after 
N jumps and P x (0) — 8 X \. The Master equation describ- 
ing the discrete time problem is given by 

P (N + 1) = Qi(l)i\(JV) + P (N) 

since the origin is absorbing 

P 1 {N + 1) = Q L (2)P 2 (N) 

P 2 (N + 1) = Q fl (l)P x (iV) + Q L (3)P 3 (N) 

P X (N + 1) - Q R (x - l)P x -x{N) + Q L {x + l)P x +i{N) 

P^iN + 1) = P L (N) + Qr(L - 2)P L _ 2 (N) 

P L (N+l) = Q R (L-l)P i _ 1 (N). (24) 

The probability to be absorbed for the first time at x = 
after N + 1 jumps (the discrete time) is 

F d[s (N+l) = Q L (l)P 1 (N). (25) 

The discrete survival probability is given by 

S dis (N) = l-P (N). (26) 

Using Eq. |22J) 

%sW = 1 - IQl(1)Pi(N - 1) + P (^V - 1)] , (27) 
and from Eq. (|2*51l 

S dis (N) = 1 - [F dig (AO + P (N - 1)] . (28) 
Using Eq. I|26|l we have 

Sdis( N )-S dis (N-l) = -F dis (N) 7 (29) 
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which simply means that the change in the survival prob- 
ability at step N is equal to minus the probability of first 
passage. Using the z transform Eq. $Z2fy of Eq. $Zj$ we 
find 



%s ( z ) = 



1 - z 



Hence from Eq. 

^CT(«) = ^{l-^dis [^(«)]} 



(31) 



Let / CT (t) be the first passage time PDF of the CTRW 
problem. As-usual 



ICT (t) = --S CT (t). 



(32) 



which is the continuous pair of Eq. (|29|) . If the random 
walker always returns to the origin, then J °° fcT(t)dt = 
1, and Eq. yields 



fcr (u) = ~uS C t(u) + 1 
and using Eq. i|3T|) 

fcr (u) = F dis Vi (u) 



(33) 



(34) 



This is the most important equation of this sub-section. 
At-least in some cases the solution of the discrete time 
first passage time problem, in z space is possible, and 
then we can transform the solution to Laplace u space 
of the seemingly more difficult case of continuous time. 
Note that our assumption that the random walk is re- 
current is valid only when the system size is finite, and 
Ql{x) > for any x besides on the boundary. 



B. First Passage Time for Unbiased Case 

We now find the first passage time distribution for the 
unbiased CTRW in Laplace space. For the unbiased ran- 
dom walk we have Ql{x) — Qr(x) = 1/2, for x ^ 0, 
x 7^ L. And as mentioned x = is the absorbing bound- 
ary condition, while L is a reflecting wall. As shown we 
may consider the first passage time for the discrete time 
random walk Eq. I|24|l and then use the transformation 
Eq. (|34|l to obtain the corresponding CTRW first pas- 
sage time. Using Eq. I|22(l the z transform of Eq. 124|) 
is 

Po («) = - Pi (z) + zPo (z) 
using the initial conditions -Pi(O) = 1, 

Pi (*)-! = \h (z) , 
for x = 2, • • • , L - 2, 



P x (z) = - P x ~i (z) + P x+1 (z) 



P L _, (z) = zP L (z) + -P L _ 2 (z) 



P~l (*) = \h-i ( z ) . 



(30) and using Eq. 



(35) 



(36) 



To solve these equations we use a recursive solution 
method [38|,|3!j. We define <f) x (z) using the relation 

P x {z) = <f) x {z)P x _ 1 {z), (37) 

and it is easy to show using Eqs. (|35l I37fl 

<f> L (z)=z/2 cj )i _ 1 (z) = (z/2)/(l~z 2 /2). (38) 
The function (j) x (z) also satisfies the recursion relation 

(z/2) 



1 - z(p x (z) /2 
which is easy to obtain from Eq. (|35|l . Let 

9x (z) 



4>x (z) 



and using Eq. I|39|l 



9x-i (z) 



h x -i (z) 



h x (z) 



9x (z) 
-f I )\h x (z) 



(39) 



(40) 



(41) 



Since we are interested only in the ratio g x (z) /h x (z) we 
may set (z) = 1 and gi(z) = z/2 using Eq. 
Eq. I|38|l gives the seeds for the iteration rule Eq. H41[) : 
hi_ 1 (z) = 1 — z 2 /2 and gi_ 1 (z) = z/2, which yield 

h L-2( Z )'9L-2( Z ) etc - Let 



h x [z)=B + (A + ) L - X + B_ (A_) L - X 



(42) 



and from h;(z) — 1 we have B + + B- = 1. A± are eigen 
values of the matrix in Eq. 1)41(1. 



1 ± vT 



A- 



Using , (z) = 1 — z 2 /2 it is easy to show 



1 - z 2 /2 - A H 



A_ - A_i 



(43) 



(44) 



and B + = 1 — _B_ . Using 
Pi(z) 



1 - z0 2 (z) /2 

and 4> 2 {z) = zh 3 (z)/2h 2 (z) and Eqs. pplli2) we find 



(45) 



4w 



z/2 



(46) 



z 2 a + J \ + 
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This equation is important since it yields the discrete first 
passage time probability with which the CTRW PDF of 
first passage time can be obtained. Using the Laplace 
transform of the waiting time PDF Eq. (|f4|l and Eqs. 
we obtain the small u behavior 



/ cx (u) ~ 1 - (2L - l)Au a + 



(47) 



To summarize Eq. (|47H yields the Laplace transform of 
the first passage times of the unbiased CTRW with re- 
flecting boundary condition on L, absorbing on the ori- 
gin, and initial location of the particle oni= 1 . 



L ) h 




C. First Passage Time for Uniform Bias 

We now find the first passage time distribution for the 
biased CTRW in Laplace space, skipping many of the 
algebraic details. Now the probability to jump left is 
Ql{x) = q and hence the probability to jump to the 
right is Qr{x) = 1 — q, for x ^ 0, x ^ L. The two 
boundary conditions are: x — is absorbing, while L 
is a reflecting wall. Like the unbiased case we treat the 
problem of the discrete time random walk and then use 
the transformation Eq. I|34|l to obtain the corresponding 
CTRW first passage time distribution. 

In this case the z transform of the master Eq. I|24f> is 

P {z) = zqP l (z) + zP (z) 
Pi (z) - 1 = zqP 2 (z) 
P x (z) = z(l-q) P x _i 0) + zqP x+1 0) 
P L _ 1 (z) = zP L (z) + z(l-q)P L _ 2 (z) 



P L (z) = z(l-q)P L _ 1 (z). (48) 

And using Eq. 

P dis (z) = zqP 1 (z) . (49) 

The solution of the biased master equation (|48fl follows 
the same procedure as for the unbiased and yields 



F(z) 



qz 



1 — q (I — q) z' 



B + \ I r 2 +B-\ L r 2 



(50) 



where 

A± (z) = 
B + + B- = 1, and 

B+{z) 



1 ± v/T= 4gz 2 (1 - q) 



1 - A- - z 2 (I - q) 
A+ - A_ 



(51) 



(52) 



FIG. 1: Boltzmann's equilibrium for an ensemble of CTRW 
particles in an harmonic potential field, and fixed tempera- 
ture. In simulations (cross) the CTRW particle with a = 
0.3, 0.5 and a = 0.8 was considered. The figure illustrates 
that for ensemble of particles, standard equilibrium is ob- 
tained, ergodicity breaking is found only when long time av- 
erages of single particle trajectories are analyzed. The scaled 
potential (dot dash curve) is the harmonic potential field, and 
the theoretical curve is Boltzmann equilibrium distribution. 
To construct the histogram we used iV = 10 6 particles, tem- 
perature T = 3, and the total observation time t = 10 6 . 



Using Eq. (|49() one can show that F^ s (z = 1) = 1, for 
any finite L and q ^ 0, namely if we wait long enough the 
particle always reaches the sticky boundary on x = 0. 

We now return to the CTRW problem. We use the 
relation Eq. I|34l) and insert in Eq. H5U|) the small u 
behavior of the Laplace transform of the waiting time 
PDF Eq. lfTl|) . In the limit of u — > we find the Laplace 
transform of the PDF of the first passage time of the 
CTRW particle 



for (u) ~ 1 - 



Au a 
2q-l 



l-2(l-g) 



L-l 



This is the main result of this section, since it will yield 
the non-ergodic properties of the biased CTRW. We see 
that for q = 1 or L = 1, Jct («) ~ 1 — Au a as expected 
since then /ct{u) = ip{u). The second term on the right 
hand side of Eq. (|53|l will diverge when q < 1/2 and 
L — > oo, as expected for an infinite system, and for a 
random walker moving against the average drift. We 
see from Eq. (I53|) . that the PDF of first passage times 
/ctW ^ t~( 1+a \ in the limit of long times, when a < 1. 
In the limit q — > 1/2 the solution for the biased case Eq. 
(|53|l . reduces to the unbiased solution Eq. I|47|l 
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V. MAIN RESULTS 
A. Non- Thermal random walks 

First consider the unbiased one dimensional CTRW 
on a lattice x — — L, • • • , L. The PDF of the fraction of 
occupation time p x = t x /t on a lattice point x, excluding 
the boundary points, is obtained using Eqs. ( Q7J 
IHH 1871 ) . The general idea of the proof is to note that 
4> x (u) = ^{u) 

i>n*(u) ~ 1- A(2L- l)u a (54) 
for u — ► and hence using Appendix A we find 

lim f{p x ) = 5 a ((2L-\)-\p x ). (55) 

Where the 8 a function was defined in Eq. Eq. l(55*)l 
does not depend on the position x of the observation 
point, reflecting the symmetry of the problem. From 
Eq. (|55Jl we see that the amplitude ratio satisfies 7l x = 
1/(2L — 1) < 1 when L > 1. This inequality means that 
we are less likely to find the particle on the particular 
lattice point x under observation (state x), if compared 
with the probability of finding the particle on any of the 
other lattice points (state nx). 



For the biased random walk, when the probability of 
jumping left is q, we consider the PDF of fraction of time 
p x on a lattice point x. Now clearly different locations 
have different distributions of the fraction of occupation 
time, reflecting the fact the the system is biased. The 
Laplace transform of the sojourn times on x is simply 

$ x (u) = $ (u) ~ 1 - Au a (56) 

and using Eq. Q19[) the sojourn times in all other states 
(nx) is 

4>nx (u) = (1 - q) fa (L -x,u) + qf L (x + L,u). (57) 

Here fn(L — x,u) is the Laplace transform of the first 
passage time PDF, for a system of size L — x+1, obtained 
in Eq. i|53|) . Similarly for (L — x,u) however now 
replace q with 1 — q in Eq. J53Jl. Using Eq. (|57|) we find 
the small u behavior of ip n x{u), and then using Eq. I)84|) 
we find 

$ nx (u)~ 1- + (5 8 ) 



Tlx = 



2q-l 



1-9 



L+x-l 



-(i-«r 



1-3 



L-x-l 



- 1 



(59) 



The latter Eqs. (|58I59|I and the results obtained in Appendix A indicate that the PDF of fraction of occupation time 
is 



with 7Z X given in Eq. I|59() . 



f(Px) = S <* (Rx,P x ) 



(60) 



As expected the PDF of the fraction of occupation time, 
for the biased CTRW, depends on the location of the 
site under consideration. As-usual if q < 1/2 the particle 
prefers to stick to the right wall. In our case this behavior 
implies that if g < 1/2 and x ~ — L (L is large) then 
Tlx — * 0, which means that the lattice point x is never 
occupied, as expected. 



B. Equilibrium— Ergodicity Breaking Relationship 

Eqs. 1(551 IdTTI ) describe the non-ergodic properties 
of the CTRW for biased and unbiased cases. We will 
now consider a relation of the problem of non-ergodicity 
with the equilibrium of the process. Consider an en- 
semble of independent random walkers performing the 
CTRW process, in the finite domain. After a long period 



of time an equilibrium will be reached, for which the den- 
sity of particles is found in a steady state profile. Such 
an equilibrium is obtained after each individual member 
of the ensemble made many jump events (one can easily 
prove that such an equilibrium is reached). We denote 
the probability of finding such a random walker on point 
x with P x q . It is straightforward to obtain P x q , though 
some care must be made when we take into consideration 
the boundary conditions of the problem. In equilibrium 

Px q = (61) 
and on the boundaries 
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peq _ 

— L — 



1-g 
9 



-L+l 



Z 



(62) 



And Z is then obtained from J2 x =-l ^x q — 1- Here Z is 
a normalization constant of the problem, not necessarily 
related to Boltzmann Gibbs statistics. 

Using the equilibrium properties of the system, after a 
short calculation of the normalization constant and some 
algebra, we find that Eqs. 1601 1591 ) may be written 
in a more elegant form 



f(Px) = S a 



peq 



l-P, 



.eq 



Px 



(63) 



0.4 



0.3 - 
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Note that P x q yields the equilibrium properties of many 
non-interacting random walkers, or the density profile 
of large number of particles. Hence the single particle 
non-ergodicity is related to statistical properties of the 
equilibrium of many particles. The fact that we find such 
a relation should be anticipated, since if we average p x 
namely consider (p x ) = J p x f (p x ) dp x we must obtain 
P x q , hence / (p x ) must be clearly related to P x q . And the 
requirement (p x ) = P* q implies that TZ X = P x q /(l~P x q ) 
as we indeed found (and similar to our discussion in Sec. 
Hill. For the unbiased case, q — 1/2 we have P x q = 1/2L, 
which leads to (JSSJ. Note that the equilibrium population 
on the boundaries x = ±L is half the value of that found 
on x ^ ±L, and hence Z = 2L even though we have 
2L + 1 lattice points. 

A possible extension of our result: we believe that if 
we consider the occupation times on M < 2L + 1 lattice 
points, Eq. © is still valid and P% = M/(2L) when 
L is large. A proof of iJBJ based on the calculation of 
the first passage time for such a case is cumbersome, if 
we consider a general configuration of M lattice points 
under observations, however we did verify this result nu- 
merically. 



C. Thermal Random Walks 



FIG. 2: The PDF of occupation times p x = t x /t where t x is 
the total time spent on lattice point x = 0, the minimum of 
the harmonic potential field, and t is the measurement time. 
Here we use a — 0.3. For an ergodic process satisfying de- 
tailed balance, the PDF f(p x ) would be delta centered around 
the value predicted by Boltzmann which is given by the ar- 
row. In a given numerical experiment, it is unlikely to obtain 
the value of p x predicted by Boltzmann, though Boltzmann 
statistics does yield the average of p x over many measure- 
ments. The PDF has a U shape indicating that events, where 
the particles hardly ever occupies x — or nearly always oc- 
cupies x = are important. To construct histograms we used 
10 6 trajectories, measurement time t — 10 6 , and temperature 
T = 3. The solid curve is the analytical formula Eqs. |S||7| 
used without any fitting parameters. 



+ Simulation 
— Analytic 




0.2 0.4 0.6 O.t 



If the CTRW particle is interacting with a thermal heat 
bath, we can relate the non-ergodicity to Boltzmann- 
Gibbs statistics. For the free particle we recall that Boltz- 
mann probability of occupying a lattice point is simply 



r 1 
P B = — 
Z 



(64) 



and as mentioned Z = 2L is the normalization condition, 
or the partition function of the problem. Hence rewriting 
Eq. l|B3|) 



f(Px)= s » 



l-P? 



iPx 



(65) 



The factor P^ /{I - ) means that with probability P x 
the particle is in state x, and with probability 1 — P x the 



FIG. 3: Same as Fig. (|5J however now a = 0.5 



particle is in state nx i.e., the rest of the system (here we 
mean probability in the ensemble sense). 

For biased CTRW when detailed balance condition Eq. 
(|T7|l holds, we find once again 



f(Px)= S » 



P l 



1-pi 3 



iPx 



(66) 



and now 



Pi 



exp(-^) 
Z 



(67) 



10 



0.03 




FIG. 4: Same as Fig. however now a = 0.8. Unlike 
Fig. J5J the PDF has a (distorted) W shape. A peak close 
to Boltzmann's value (the arrow) is an indication that as a is 
increased the ergodic phase is approached 



where x is the lattice site under observation, V(x) — 
—Tax is the potential field, and a is the lattice spacing. 
Here the partition function is 



Z = 



T 



i \-L + l /, \L-1 



2q - 1 



(68) 



which is easily verified once proper reflecting boundary 
conditions are applied, and using Eq. (|17|) . 



D. Numerical Demonstration 

In previous sections we considered the cases of biased 
and unbiased CTRWs. We see however that our results 
may be more general, and valid also for random walks 
in a general deterministic external field. We decided to 
check this issue using the example of a random walk 
in an Harmonic trap. For that aim we used numerical 
simulations, since calculations of the first passage time 
are cumbersome. The problem of anomalous diffusion in 
Harmonic potential was considered in the context of frac- 
tional Fokker-Planck equations ^(j an d m single particle 
experiments |34j . Anomalous diffusion in harmonic field 
was also investigated using fractional Langevin equations 
|4ll . It would be interesting to test if such stochastic 
equations yield an ergodic behavior. 

The potential field we choose is V(x) = Kx 2 , and K = 

1. We used: (i) the condition of detailed balance Eq. 

(|18fl . and (ii) at bottom of the well, point x = 0, we 

used the symmetry of the potential and choose Ql(0) = 

Qr(0) — 1/2- These two conditions yield Ql(x). In 

simulations we generate random waiting times, according 

to the normalized power law waiting time PDF ip(r) — 
aT -(i+ a )^ for T>1 

We first checked that Boltzmann equilibrium is reached 
for an ensemble of particles. In these simulations we build 
histograms of the position of N = 10 6 particles, after 



FIG. 5: The PDF of fraction of occupation time, for a = 0.8 
when temperature T is varied. Similar to Figs. 0-Elwe con- 
sider occupation times on point x = of the Harmonic poten- 
tial. When temperature is such that probability of occupying 
x = is equal 1/2 the PDF is symmetric, for our parameters 
this temperate is T s — 1.3. For T « T a the particle tends 
to be located all the time on x = 0, hence for T = 0.8 < T a 
(the dot-dashed curve) the PDF has more weight on values 
p x > 1/2. For T >> T s the particle is never found on x — 0, 
hence for T = 3 > T s (dotted curve) the PDF has more weight 
on values of p x < 1/2. 



each particle evolves for a time t — 10 6 . In Fig. \I\we find 
good agreement between our simulations and Boltzmann 
statistics when many particles are considered. The Fig. 
illustrates that an observer of a large number of particles 
cannot detect ergodicity breaking, and the single particle 
limit is essential for our discussion. 

We then consider one trajectory at a time. We obtain 
from the simulations, the total time spent by the particle 
on lattice point x = 0, namely at the minimum of the 
potential. This time is t x and the fraction of occupa- 
tion time p x — t x /t. In the ergodic phase and long time 
limit p x will approach the value predicted by Boltzmann 
statistics. While in the non-ergodic phase we test if our 
prediction Eq. (|b17|) hold. In Figs. |3 El HI we consider 
three values of a, a = 0.3, 0.5, 0.8 and fix the tempera- 
ture T. All figures show an excellent agreement between 
our theoretical predictions Eqs. (I(il7l) and numerical sim- 
ulations. It is more important however to understand the 
meaning of the figures. 

For small < a << 1 we expect that the particle will 
get stuck on one lattice point during a very long period, 
which is of the order of the measurement time t. This 
trapping point, can be either the point of observation 
(e.g. x = in our simulations) or some other lattice 
point. In these cases we expect to find p x ~ 1 or p x ~ 
respectively. Hence the PDF of p x has a U shape. This 
case exhibits large deviations from ergodic behavior, in 
the sense that we have a very small probability for finding 
the occupation fraction close to the value predicted based 
on Boltzmann's ergodic theory. As shown in Figs. 
such U shape behavior is found for the cases a = 0.3 
and a = 0.5. We also plotted the prediction made using 
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the ergodic assumption (the arrows in the Figures) to 
demonstrate the fact that a measurement is not likely to 
yield the average which is located on P x . 

When we increase a we anticipate a more ergodic be- 
havior, in particular in the limit a — > 1. An ergodic 
behavior means that the PDF of the occupation fraction 
p x is centered on the Boltzmann's probability. In Fig. 
21 where a — 0.8 we start seeing a peak in the PDF of 
p x centered in the vicinity of the ensemble average value. 
Note however that the PDF f(p x ) still attains its max- 
ima on p x = and p x — 1 . Hence we find a weaker non- 
ergodic behavior if compared with the cases a = 0.3, 0.5, 
and a W shape of the PDF. 

In Fig. [5] we consider the example of a particle in an 
Harmonic field, we fix a = 0.8 and vary temperature, us- 
ing Eqs. I|6I7|) . The observation point remains x = 0. At 
temperature T ~ 1.3 (solid line) we see that the PDF of 
p x is symmetric. This happens when P x = 1/2, namely 
for a case that there is probability half of occupying the 
observation point, and probability half to be out of this 
point. When the temperature is very low, we expect to 
find the particle, in the ground state, namely on x = 0. 
Hence the PDF of p x is tilted towards p x ~ 1 when tem- 
perature is lowered (see Fig|S]wlien T = 0.8 ). In contrast 
when the temperature is high, we expect the probability 
of occupying the observation point x — to be reduced 
(as-usual entropy wins at high temperature) . And indeed 
we observe that when T — 3 the PDF of p x is more tilted 
towards the left namely to p x ~ 0. 



E. Validity of main Eqs. (15171) 

Our numerical work as well as our analytical solutions 
for the bias and non-bias CTRW show the validity of 
Eqs. (|6I7|I . What happens for more general type of po- 
tential fields? Can we claim, that Eq. @) has a wider 
applicability? Consider the CTRW with potential profile, 
{...E x ...}, with the dynamics satisfying detailed balance 
condition. We claim, but have no rigorous proof, that if 
for ip{ T ) with finite moments, the system is ergodic, then 
for the same energy profile but when the waiting time 
has a long tail, Eqs. H6I7[) hold. Our reasoning is that 
we can think of a as a control parameter, which we can 
vary between < a < 1 . And since for the case a = 1 we 
have TZ X = P x /{I - P x ) also for < a < 1 this relation 
must hold (since Pg does not depend on a). Further the 
transformation s(t) — > As" in the small s behavior of 
the waiting time seems to indicate that the behavior we 
found has a general validity. 

A way to understand the ergodicity breaking laws, Eqs. 
(|6I7() . is to consider the number of times n x the particle 
visits lattice point x, during a long measurement time. In 
that case the particle visits x many times, and we assume 
that the fraction of number of visits satisfies 



where n is the total number of jumps made by the parti- 
cle. If the first moment of the waiting time distribution 
is finite, we have t x /t — n x /n, since the average time 
spent on x is n x times the mean waiting time. When 
the average waiting time is infinite a < 1, one can show 
that the PDF of t x /t is given by Eq. © if condition Eq. 
(|69|l holds. Eq. (|69|) should be tested in more detail, for 
example using numerical simulations. 



VI. SUMMARY AND DISCUSSION 

We obtained the non ergodic properties of biased and 
unbiased continuous time random walks. In particular 
the distribution of the occupation fraction p x was found. 
Our results are valid both for thermal and non-thermal 
cases. In both cases the non-ergodicity is described using 
the 5 a (lZ x ,p x ) PDF Eq. JBJ). Where a is the anomalous 
diffusion exponent (a; 2 ) ~ t a . Both for thermal and non- 
thermal random walks the parameter 1Z X is related to the 
ensemble averaged equilibrium properties of the system, 
Eqs. (|63l 1661) respectively. If the system is in vicin- 
ity of thermal equilibrium, the equilibrium of the sys- 
tem is the Boltzmann-Gibbs equilibrium, in the ensem- 
ble sense. Such behavior is found when detailed balance 
conditions are applied. In this case the characterization 
of the non-ergodic properties of the occupation times is 
related to the partition function, and temperature. The 
non-ergodicity manifests itself when time average of sin- 
gle particle observables is considered. In particular the 
occupation time, in a given energy state, or on a partic- 
ular lattice point. Hence the non-ergodicity might reveal 
itself in single particle experiments. 

Models and systems describing anomalous diffusion are 
wide spread. In most cases ensemble average properties 
of such processes are investigated, both in theory and in 
experiment. For single particle experiments, where the 
problem of ensemble averaging is removed, we may ei- 
ther: i) reconstruct the ensemble averages, by repeating 
the single molecule experiment many times, or ii) inves- 
tigate the ergodic properties of the system, by consider- 
ing fraction of occupation time in a particular state, and 
obtaining its distribution. It is the second type of mea- 
surement which is considered here, which yields insight 
into single particle properties which differ from the stan- 
dard ensemble measurement, provided that a non-ergodic 
phase is investigated. And while the theory of anomalous 
diffusion processes is now vast, the non-ergodic properties 
of such processes are still not well understood. Investi- 
gation of this topic beyond the CTRW approach is left 
for future work. 
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FIG. 6: A schematic diagram of the two state process. 



VII. APPENDIX A: THE S a (TZ x ,p x ) FUNCTION 

In this Appendix we re-derive the limit theorem Eq . 
© . While this goal was accomplished long time ago [3g , 
we believe that it is worth while re-deriving this result 
using methods similar to those used today in statistical 
Physics community [26|. We also derive an exact distri- 
bution for the occupation fraction of a two state process 
in Laplace space Eq. I|81l I82|) . 

Consider a system evolving between two states, + and 
— , corresponding to states x and nx respectively. Let 
6 {t) = 1 when the system is in state + otherwise 6(t) = 
and the system is in state — . A schematic diagram of 9(t) 
is shown in Fig. Let {i,} denote dots on the time axis 
on which transition events between state + to state — or 
vice versa is made. Let {t{\ be sojourn times either in 
state + or state — . If the process starts with state +, then 
Tj is a + state if i is odd. We also denote the total number 
of jumps, in the measurement time interval (0, t) with 
n. We assume that the sojourn times are independent 
identically distributed random variables. The PDF of 
sojourn times is ip+{r) and ip-i T ) f° r states + and — 
respectively. Such a simple process is called a two state 
renewal process. As-usual it is convenient to analyze such 
a stochastic process using the Laplace transforms 



4>± 0) = 



r ip± (r) dr. 



(70) 



While we will consider general properties of the 
stochastic process, we will eventually focus on two main 
cases. First consider the case where all moments of ip±(r) 
are finite, e.g. exponential PDFs belong to this category. 
Then the following small s expansion holds 



4>±(s) ~ 1 - s(t±) + 



(71) 



Here (t±) are the mean sojourn times in states ±. A 
second generic case is: 



ip±(s) ~ 1 - s a A± 



(72) 



with < a < 1. For example the one sided Levy PDF 
tp±(s) = exp(— A±s a ) belongs to this class. In the time 
domain these PDFs behave like 



A A 



\T{-a)\t l + a 



(73) 



when t is large, namely for this family of PDFs the aver- 
age waiting times in both states diverges. 

Let t+ be the total time spent in state +, within the 
time period (0, t). Then the occupation fraction in state 



P+ = T 



t t 



(74) 



We now consider statistical properties of t+ focusing on 
the scaling limit t — ► oo. Let / t + „ (t+) be the PDF of t+ 
conditioned that n renewal (i.e. jumps) events occurred 
in the time interval (0, t), and that the start of the process 
is in state +. 

Consider the case n odd, n = 2k + 1 with k = 0, 1, • • •. 
Then since we start with state +, t+ = X^Liodd Ti 
where the summation is only over odd Vs. Also we have 
t n < t < t n+ i, where i„+i is a renewal event which occurs 
after end of measurement (see Fig. EJ). Hence 



fun (*+) = (S\t+- J2 T A 1 (*» ^ 1 ^ 

i=i,odd 



where 5 (x) is the Dirac delta function, and 



(75) 



, , > J 1 If condition in parenthesis is true 

1 {t n < t < t n+1 ) - | Q otherwige 

(76) 

In Eq. (|75l) (• • •) denotes an average over the stochastic 
process, soon to be specified in more detail. Later we will 
consider the case n even, and then sum over n to obtain 
the PDF of t + . 

Let s be the Laplace pair of t, and u of t+. It is con- 
venient to consider the double Laplace transform /+ s (u) 
of /+*(*+) 



fn,s («) = / ' 
10 



e- st f+ (t+)dt + dt = 



oo poo 



JO 



t S\t + - n\I(tn<t< t n+1 )dt A 

i=i.odd 



e "^,=i,odd r ') ; (77) 



where we made use of Eq. I|75|l . Now we may consider 
the average (• • •), using t n = £" =10( jd n + Eleven 7 "* 
and similarly for i n +i. Recalling that Ti with odd (even) 
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i are + (— ) states respectively we find after averaging must define t„ = t — t n as the time difference between 



over the {rj}s 



fls {U)=i,% +I {s + U)^_ ( S ) 



(78) 



where n — 2k + 1. For even n such that n = 2fc, fc = 
1, 2, • • • we obtain 



?4- i \ 7k , \ Ik / x 1 — V"+ ( S + U ) 



s + w 



(79) 



end of the measurement and last jump in the sequence 
(see Fig. EJ). We are ready to obtain the double Laplace 
transform of / t + (<+), i.e., the PDF of t + when the process 
starts in + state, 



00 



(80) 



Using Eqs. H78I79I8()|1 we obtain the exact solution to the 
Note that for even n the last + interval falls on t, and we problem in Laplace s, u space 



i-i>+(s + u) ? (s) 

■ h ?/>+ (s + u) 

S + It s 



1 



1 - "0+ ( s + u ) ty- ( s ) 



(81) 



It is easy to check the normalization condition f^~(u — 0) = 1/s provided of-course that ip±(r) are normalized PDFs. 
In a similar way one can show that if we start the process in state — the double Laplace transform of the PDF of t + 
denoted with (t+) is 



(«) 



V- (a) 



1 - 1p + (S + It) 1 - 7/L (s) 



s + u 



1 — ip + (s + u) ip- (s) 



(82) 



Eqs. (|81l I82|l yield in principle the exact expression for 
the occupation fraction, which might be useful in deter- 
mining the pre-asymptotic behavior, for example using 
numerical inverse Laplace transform. 

For the generic case Eq. 11721 . in the limit of s — > 
and u — ► their ratio remaining arbitrary, Eqs. I|81I82() 
yield 



ft («) 



K(s + u) L 



n(s + u) a + s c 



with 



K = 



A J 



(83) 



(84) 



The amplitude ratio 1Z determines the degree of symme- 
try in the problem. Note that in this scaling limit the 
initial state of the process, i.e. process being in state + 
or — at initial time, is not important. 

The small (s, u) limit considered in Eq. (|83|l corre- 
sponds to large measurement time t, and occupation time 
t + limit. We invert Eq. iJSSJl using a method given in 26] . 



The method states that if in the limit s, it 
Laplace transform behaves like 



a double 



/. (u) = -9 (-) (85) 
then the PDF of the scaled variable p + = t + jt is in the 



long time t limit 



f (p+) = limlmg 



1 



x + ie 



(86) 



Using Eq. (|83(l we find the PDF of the fraction of occu- 
pation time p + = t+ jt 

f (p+) = S a (K,p + ) = 



sm Tree 



m- l 0--p + ) a ~ 1 



n n 2 (i-p. 



,2a 



■p 2 ? + 2TZ (l - p + ) a p% cos it a 

(87) 

The PDF is normalized according to J Q f(p + )dp + = 1, it 
is valid only in the long time t limit and is independent 
of it. In this sense an equilibrium is obtained. In par- 
ticular when A + = A- and a = 1/2 we find the arcsin 
distribution. It is easy to show that the average 



(p + ) = 



A, 



t A + 
In the limit a — > 1 we obtain 



A. 



(89) 



where (t)± are the average waiting times when the wait- 
ing time PDFs have finite moments. We identify this 
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behavior with an ergodic behavior, since according to 
Eq. (I74|) p + is a time average of 8(t), which is equal to 
the ensemble average value when moments of ip± (r) are 
finite. 

From Eq. (|87|l we see that the PDF of p + is not 
narrowly centered on the ensemble average value when 
a < 1. Hence the case a < 1 is called the non-ergodic 
phase. Indeed any measurement of p + is unlikely to yield 
the average, and / — * oo when p + — > or p + — ► 1. 
The latter behavior corresponding to systems in state — 



or state + during nearly all the measurement time, re- 
spectively. The reason for non-ergodic behavior is that 
in a measurement time interval (0, t) we expect to ob- 
tain a few sojourn times of the order of t, these can be 
either + times or — times or both. In each measurement 
we make these large sojourn times will be different than 
those found in a second measurement. Hence the time 
average of 0(t), the occupation fraction p, will remain 
random even in the long time limit. 
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